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Abstract 

The irreducible representations of the Lie algebra su(3) describe rotational 
bands in the context of the nuclear shell and interacting boson models. 
The density matrices associated with su(3) provide an alternative theoretical 
framework for obtaining these bands. The su(3) density matrix formulation 
is mathematically simpler than representation theory, yet it yields similar 
results. Bands are solutions to a system of polynomial equations defined 
by the quadratic and cubic su(3) Casimirs. Analytic solutions are found in 
many physically important cases including rotation about principal axes and 
spheroids. Numerical solutions are reported in other cases including tilted ro- 
tors. The physics of su(3) rotational bands is more transparent in the density 
formalism than in representation theory. 
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I. INTRODUCTION 



The physics of many-particle systems is dominated in many cases by a relatively small 
number of degrees of freedom. The quintessence of the physics is clarified by models which 
ignore minor effects and focus on the principal dynamical factors. When the set of relevant 
observables closes under commutation to form a Lie algebra, it is called a spectrum gener- 
ating or dynamical symmetry algebra [?,?,?], and a group theoretical model is suggested as 
the natural explanatory framework. 

If g is a spectrum generating Lie algebra of hermitian operators, then its irreducible 
unitary representations (irreps) define group theoretical models. The decomposition of the 
reducible representation of q on Fock space into its irreducible subspaces provides the micro- 
scopic interpretation of the models. Although an irrep of q simplifies the original Fock space 
problem, its dimension may be still too large to allow an easy analysis of the physics, e.g., 
noncompact algebras like sp(3, M) of the symplectic collective model [?,?,?] or u(6, 6) of the 
extended interacting boson model [?] have infinite-dimensional irreducible representations. 
Intractable models also may arise when the rotation group algebra so (3) is not canonically 
embedded in q and the angular momentum decomposition is not multiplicity free, e.g., the 
compact symplectic group sp(2j + 1) which describes seniority in a single j-shell [?,?], or 
the s<ig-boson model based on u(15) [?,?]. 

The relationships among the shell model, Hartree-Fock, and the group U (n) of unitary 
transformations in the n-dimensional single-particle space suggest a solution [?,?,?]. 

The shell model may be regarded as a group theoretical model in which the spectrum 
generating algebra is the set of all one-body hermitian operators u(n). Throughout this 
article, Lie algebras are denoted by jj, e.g., u(n), su(3), so (3), while the corresponding Lie 
group is written as G, e.g., U(n), SU(3), SO(3). Each Lie group is the exponentiation of its 
Lie algebra, e.g., U(n) = exp(u(n)). For k identical valence fermions, the shell model space is 
the totally antisymmetric irreducible representation of U(n) with dimension n\/[k\(n — &)!]. 
Already for medium mass nuclei with active valence neutrons and protons, the factorially- 
growing dimension of the antisymmetric irrep is astronomical. Moreover the rotation group 
is not canonically embedded in U(n). 

A way around this difficulty is self-consistent mean field theory. In Hartree-Fock, the 
mean field Hamiltonian must be diagonalized in the single-particle space of dimension n - 
no matter how many particles are in the valence space. The set of admissible states in the 
Hartree-Fock approximation is the set of Slater determinants, or, equivalently, the set of 
idempotent hermitian density matrices, p 2 = p, with trace equal to the number k of valence 
nucleons, tr p = k. The density matrices are defined by the expectation of the one-body 
hermitian operators in u(n). Mathematically, the set of density matrices spans the dual 
space of the algebra u(n). The set of admissible idempotent densities is not a vector space, 
but a surface of dimension 2k(n — k) contained within the vector space of all hermitian 
density matrices [?]. 

The surface of idempotent densities is a level surface of the Casimirs. There are n 
Casimirs for the unitary algebra u(n) given by the trace of powers of the density, C r = tr p r 
for r = 1, ... ,n. On the idempotent density surface, the Casimirs are evidently constant 
C r = k for all r. Conversely, the level surface in the dual space of the Casimirs C r = k for 
integral k G [0, n\ consists of the Hartree-Fock densities. 
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The idempotent densities transform among themselves by the elements of the unitary 
group: a unitary matrix g transforms a density p into gpg~ x . U (n) is a transformation group 
on each level surface of the Casimirs. 

Thus there is a common algebraic structure for both the shell model and Hartree-Fock. 
The distinctions arise from the way the algebra determines the model states. For the shell 
model, the quantum states span an irreducible unitary representation of the group U (n). For 
Hartree-Fock, the role of the unitary group and its Lie algebra has three essential aspects: 
(1) the densities in mean field theory are elements of the dual space of the Lie algebra, (2) 
the allowed Hartree-Fock densities are a level surface for the Casimirs, and (3) the unitary 
group transforms the allowed densities among themselves. These are the three ingredients 
for a mean field theory that may be implemented for any Lie algebra. 

The aim of this article is to construct a density matrix or mean field theory for the test 
case of the 577(3) model that underlies the algebraic theory of nuclear rotational motion. 
SU (3) is a paradigm algebraic theory since it is the simplest nontrivial dynamical symmetry 
in nuclear structure physics. The relationship between SU (3) irreducible representations and 
their mean field approximations is investigated. It is shown that there is a close correspon- 
dence between the results from these two models. But the 577(3) density matrix theory is 
mathematically simpler and provides a clear physical interpretation of the rotational bands 
found in the irreps. 



II. 611(3) DENSITY MATRICES 

A. Algebra definition 

Let {x a j,p a j) denote the Cartesian components of the dimensionless position and mo- 
mentum hermitian operators of particle a in a system of k particles. They obey the canonical 
commutation relation [x a j,p/3k] = i5 a pbjk- The traceless Elliott quadrupole operator [?] 

Qjk 2 (^oj^ak PajPak 3 $jk (%afi%afi Pa^Pafi)^) (^-1) 

and the vector angular momentum operator 

Ljf- %ajPak ^akPaj 
Li — \ SijkLjk (2-2) 

(summation over repeated indices) generate an eight-dimensional real Lie algebra of one- 
body hermitian operators, 



[ Lj -i Lfc\ i £jkmL m 



[Qfk > 2-] = i 1 (s kr e jsa L a + 5 ks e jra L a ( 2 - 3 ) 

-\-8j r £ k sa L a -\- 5j s Ek ra LaJ , 



that is isomorphic to the algebra of hermitian traceless matrices 
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su(3) = {Z e M 3 (C) | = Z, tiZ = 0} 



(2.4) 



If X and y are real 3x3 matrices, Z = Y+iX G su(3) if and only if X T = -X, Y T = Y, and 
tr y = 0. The isomorphism a between the algebra of matrices and the algebra of hermitian 
operators is given explicitly by: 

o{Z) = Y jk Qf k - \ X jk L jk (2.5) 

for Z = Y + iX G su(3). Note that [a(Z), a(W)} = i a([Z, W\) for Z, W G su(3). 

Alternatively, the algebra of hermitian operators may be defined as the su(3) dynamical 
symmetry algebra of the interacting boson model [?,?,?]• If S and denote the boson 
creation operators, then the generators are the angular momentum and quadrupole operator, 

$? = fy + 8%^[<fx^\ (2.6) 

where =F corresponds to particle and hole bosons, respectively. This IBM algebra of hermitian 
operators is also isomorphic to the algebra of matrices su(3). 

Mathematically, all isomorphic copies of an algebra are indistinguishable. The su(3) 
algebra of matrices is more convenient to use for calculations than either the Elliott or IBM 
operator algebras. 



B. Dual space 

Consider an algebra of hermitian operators acting on the Hilbert space Ti which is iso- 
morphic to su(3), e.g., the Elliott or IBM representations. For each normalized state vector 
^ G H the expectations of the operators 

q jk = <tf | Qf) | *) 

l jk = (*\L jk \iit). (2.7) 

define a real symmetric traceless matrix q and a real antisymmetric matrix /. The "density" 
matrix corresponding to ^ is defined as the hermitian traceless matrix p — q — ^il. In terms 
of it, the expectation of a general element of the operator algebra is 

(p,Z) = tr(pZ) = <tt|<7(Z)|tt), (2.8) 

for Z = Y + iX e su(3). 

In fact each traceless hermitian density matrix p defines a real-valued linear functional 
on the matrix Lie algebra su(3), viz., (p , Z) = tr(pZ) for all Z G su(3). The set of all 
such linear functionals is called the dual space of su(3) and is denoted by su(3)*. In Dirac 
quantum mechanics the dual space is the space of "bras." The mapping from the Hilbert 
space to the dual space is called the moment map M : TL — > su(3)* where the density 
corresponding to the vector \l> is p = M(\&) [?,?,?]. 

The density retains only part of the entire information about the system that the wave 
function carries, but a very important part - the expectations of the su(3) observables. 
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The dimension of the dual space is the same as the dimension of the su(3) algebra and is 
a significant simplification of the quantum problem in the Hilbert space 7i. It reduces all 
the degrees of freedom incorporated in the wave function to just those most relevant to the 
physics of su(3) rotational states. 

When ty HW is a highest weight vector for an irreducible representation of su(3), the 
corresponding density is a diagonal matrix. To see this, express the Elliott generators in 
terms of the one-body operators 

Cjk = \ (a^a^ + a ak al^j , (2.9) 

where a ak , a) ak are the harmonic oscillator bosons, a ak = (x ak + ip ak )/V2, a\ k = {x ak — 
iPak)/V%- The Elliott su(3) generators are given by 

Qfk = 2 (^jk + C kj - -5 jk H^j 

L jk = -i (C jk -C kj y (2.10) 

where H is the harmonic oscillator Hamiltonian, H = Cjj. By definition, the highest 
weight state is annihilated by the raising operators and is an eigenvector of the su(3) Cartan 
subalgebra, 

Cjk^Hw = 0, when j < k 
(C3 3 -C 11 )V HW = \y HW (2.11) 
(C U -C 22 )^ HW = /i^W, 

when the weights A,/i are nonnegative integers. Therefore the density of the highest weight 
state is 

l /-X + fJi \ 

P(x,) = ~ -A-2/i . (2.12) 
6 \ 2X + /I J 

A similar argument for the IBM su(3) algebra yields the same diagonal density matrix for 
the highest weight vector. Indeed the derivation is independent of the specific realization of 
the su(3) operator algebra. 

C. SU(3) group transformation 

The group SU(3) consists of the complex 3x3 unitary matrices with unit determinant. 
By exponentiation, a representation a of the Lie algebra su(3) extends to a representation, 
also denoted by a, of the group SU (3). Even when the Lie algebra representation is known, it 
is very difficult to determine explicitly the group representation. However the corresponding 
group transformation of the densities is simple. 

Suppose is a normalized vector in the Hilbert space that carries the unitary represen- 
tation a of SU(3). Let p = M(\l>) denote its corresponding density in the dual space. The 
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group SU(3) transforms ^ into <j(g)^ while the density is transformed into M(a(g)^f). For 
the transformed density we have 

(M(a(g)*),Z) = (a(g)V | a{Z) \ a(g)V) 

= (* | a(g)- 1 a(Z)a(g) | *) 

= (^|a(( 7 - 1 ^)|^)=tr(p( ? - 1 ^) 

= (gP9~\Z). (2.13) 

Hence the density p transforms into gpg~ l , a product of three matrices. The group trans- 
formation in the dual space is called the coadjoint action and it is denoted by Ad*p = gpg^ 1 

[?]■ 



D. Casimir Invariants 

The Casimir invariants, or Casimirs, are polynomials in the algebra generators that 
commute with all Lie algebra elements. su(3) has two independent Casimirs of quadratic 
and cubic orders, 

C r (p) = trp r , for r = 2,3. (2.14) 

These are functions on the dual space that are invariant with respect to the coadjoint 
transformation, C r (Ad*p) = C r (p) 

When p — q — the invariant functions are 

C 2 (p) = tiq 2 -\til 2 

C 3 (p) =trg 3 -|tr(g/ 2 ). (2.15) 

In particular the values of the Casimirs at the diagonal density corresponding to a highest 
weight vector are evaluated to be 

C 2 (P(a m )) = | (A 2 + Ap + p 2 ) 

C-Mxn)) = |(2A 3 + 3X 2 p - 3Xp 2 - 2p 3 ). (2.16) 

In quantum mechanics, observables are hermitian linear operators, while in density ma- 
trix theory, observables are real-valued functions on the dual space. The expectations of 
the Casimir operators with respect to a highest weight state differ from the density matrix 
functions C r (p) by terms of lower degree. The ultimate reason for the discrepancy is that 
quantum fluctuations are not included in the density matrix theory. A mathematically rigor- 
ous presentation of the relationship between polynomial operators in the enveloping algebra 
and functions of the density matrices is given in the Appendix. 

It is important to maintain consistency within the density matrix theory and not replace 
the values of the Casimir functions, Eq. (|2.16|) , by their quantum expectations in hopes of 
an improved theoretical description. Similarly the square of the total angular momentum is 
I 2 in the density matrix theory, not 1(1 + 1). 
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E. Admissible densities 



In the mean field approximation the admissible densities are restricted to those that 
lie on a level surface of the Casimirs, i.e., a surface on which the two Casimir functions are 
constant. The 577(3) group transformations Ad* leave each level surface invariant. Since any 
hermitian matrix p G su(3)* may be diagonalized by some unitary transformation g G SU (3), 
each level surface contains a traceless diagonal matrix. Because eigenvalues are unique, each 
level surface contains a unique diagonal matrix up to ordering of the real eigenvalues, which 
may be parameterized by nonnegative A and \i as in Eq. (|2.12|) . Note that a diagonal density 
only corresponds to a highest weight vector when A and \i are also integers. 

In the typical case, a level surface is six dimensional, because there are two functionally 
independent conditions imposed in the eight dimensional dual space. In the special cases of 
p or A equal zero, the level surface of admissible densities is four dimensional. 



III. 6H(3) ROTOR STATES 

The rotation group 50 (3) is a subgroup of the special unitary group 577(3). A density 
p = q — \ il in su(3)* is transformed by a rotation R G 50(3) into the density Ad* R p = 
RpR T = RqR T — \ i RIR T . Since any real symmetric matrix can be diagonalized by a rotation 
matrix, there is a R G 50(3) such that the rotated quadrupole moment is diagonal, 

q = RqR T = diag (ft, g 2) fts)- (3.1) 

The eigenvalues are unique, up to their order, which we fix to be q% > q\ > q<i- From a 
geometrical viewpoint, R rotates the laboratory frame into the body-fixed frame in which, by 
definition, the system's quadrupole moment q is diagonal . At the same time the laboratory 
angular momentum I is transformed to I = Rl R T , which is the system's angular momentum 
projected onto the body-fixed principal axes. The matrix I is antisymmetric, but otherwise 
arbitrary. In general, the angular momentum vector is not aligned with a principal axis. 

The diagonal entries of q define the (/3, 7) deformation parameters in the body-fixed 
frame [?] 

q k = /3cos( 7 -fcf), fc = l,2,3. (3.2) 

The chosen ordering for the eigenvalues q k corresponds to (3 > and 7 G [0, 7r/3]. The trace 
of any power of q is a rotational scalar; the quadratic and cubic scalars simplify to 

tr(g 2 ) = |/3 2 (3.3) 
tr(g 3 ) = §/3 3 cos(3 7 ), (3.4) 

which are model- independent measures of deformation [?]. 

The angular momentum is a pseudovector. The vector components of the angular mo- 
mentum are given by Zj = | Ij^ in the laboratory frame and by Jj = | Ij^ in the 
body-fixed frame. The rotation of the vector angular momentum I = Rl is equivalent to 
the matrix transformation I = RIR T . A principal axis rotation requires that two of the 
three components of I are zero. A tilted rotation in a principal plane requires that one 
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component of I is zero. But, in general, all three components of the angular momentum in 
the principal axis frame are nonzero. The rotational scalars of Eq. ( |2.15| ) that are quadratic 



in the angular momentum matrix may be expressed in terms of the vector components 

tr I 2 = —2 Iklk 
tr(ql 2 ) = q k I 2 k . (3.5) 



A. Angular momenta and deformations 

The range of possible angular momenta and deformations is restricted because the ad- 
missible densities lie on a level surface of the Casimirs. Since any admissible density may 
be rotated to the principal axis frame, it is sufficient to solve for the admissible body-fixed 
densities. Such densities with total angular momentum I are simultaneous solutions to the 
algebraic system: 



qi + <?2 + <?3 = (3.6) 

E <g + ^ = ° 2 ( 3 - 7 ) 

k 

k k 

II + l\ + l\ = I 2 . (3.9) 



This is the fundamental set of algebraic equations for su(3) density matrix theory. It is 
an underdetermined system of four equations for six unknowns (q±, q 2 , q%, I-i, L,L). The 
fundamental system imposes the constraint that a density is admissible and has total angular 
momentum /. 

The system (|3.6])-(pH|) determines the change in the shape of a rotating body as the an- 
gular momentum increases. Since it is a system of four equations for six unknowns (qj~, 
analytic solutions are given uniquely only when additional assumptions are imposed. Sev- 
eral important special solutions can be derived including rotation about principal axes and 
spheroidal nuclei. 

Expressed in terms of the (j3, 7) collective coordinates and spherical coordinates for the 
angular momentum in the body-fixed frame, I\ = Jcos^sin^, J 2 = Jsin^sin^, J 3 = / cos 9, 
G [0,7r], 9 e [0,tt/2], the system ([H|-([n|) is equivalent to 

3P 2 + I 2 = 2C 2 (3.10) 

3/3 3 cos(37) + 3(3I 2 A((f),e^) = 4C 3 , (3.11) 

where 

A{<t>, 9, 7) = \ ((2V3 sin 2 - y/3) sin #sin7 

+ (3sin 2 #-2)cos7) . (3.12) 

The deformation (3 is a unique function of the angular momentum /. This is a kine- 
matical property of the su(3) model which is due to the shell model compactification of the 



S 



rotational algebra, i.e., the replacement of the exact quadrupole operator by the in-shell 
Elliott expression. The energy is irrelevant to the j3{I) functional relationship. In contrast 
the triaxiality 7 and the direction of the angular momentum / are not determined uniquely 
by the fundamental system of equations. Additional assumptions are required to derive so- 
lutions, either kinematical (rotation about principal axis or spheroidal shape) or dynamical 
(energy functional E[p\). 

To these four equations, two more must be added to determine the quadrupole moments 
and angular momentum. In the following two sections, the critical points of an energy 
functional are used to supply the missing equations. But, in this section, many simple 
kinematical solutions are derived that are independent of the energy functional. Principal 
axis rotations require that two body-fixed components of the angular momentum vanish. 
Spheroidal solutions are obtained when two moments are equal. First, however, the range 
of the angular momentum that is compatible with the fundamental set Eqs. ( p.6|) -( ^9l) is 
determined. 

The function A(<p, #,7) ranges from -1 to +1. Within the chosen interval for the angle 
7, 7 G [0, 7r/3], sin 7 > and cos 7 > 0. The maximum values for sin 2 </> and sin 2 ^ are 1, so 
^4max = cos(7 — |) = 1 when 7 = ir/3. Similar argument leads to A min = —1. 

Solution III.l The minimum angular momentum is I = and the principal axis 
quadrupole moments of a nonrotating body are: 

2A + u -A + u -A-2/i 
13 = — 7, — >Qi = 7, > 92 = o • ( 3 - 13 ) 



Note that any permutation of the three axis lengths are solutions to the system 
when 1 = 0. 

The system ( |3.1U| , 1| ) has a solution for the intrinsic angular momentum / at both ends 
of the allowed intervals of values for (f), 9, and 7: I = A + fi when 7 = 0, 9 = 0, <fi G [0, ir) 
and A < \i\ I = A + \i when 7 = 7r/3, 9 = it/ 2, <fi = tt/2 and A > /i. 

Proposition 1 The maximum allowed angular momentum is the same as the upper bound 
found in representation theory, I = A + fi. 

Proof: Reductio ad absurdum. Suppose the angular momentum can have values bigger than 
A + fi. Consider the behavior of (3(1) and A(I) in ( |3.10| , |3.11| ) as the angular momentum 
increases. 

For the case of oblate spheroids 7 = 7r/3 and rotation around the symmetry axis (cf) = 
7r/2, 9 = 7r/2) the trigonometric function A = A max = 1. For A > fi the angular momentum 
is / = A + /i, and the deformation is 3/3 = A — fi. Differentiating Eq. ( |3.10|) with respect to 
the angular momentum / and evaluating the result for / = A + fi shows that as the angular 
momentum increases, the deformation parameter j3 decreases: 



dp 
dl 



A + n 
A — /U 



< 



(3.14) 



Differentiation of Eq. ([3.11 ) with respect to I and evaluating the result for I = X + fi leads 
to: 
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12 (A + /x) v _ (A-/x) 2 d(cos(3 7 )) 



(A-//) 5 
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dJ 



7=7r/3 



+( d7 

d 7 



Assuming 7(1) is a different iable function, — 7^ 00, and 

dJ 



d(cos(37)) 



dJ 



7=7r/3 



-3sin(37)^y 



0. 



7=7r/3 



Therefore 



12 7T ± 4V>0 
(A - 



(3.15) 



(3.16) 



(3.17) 



or the function A increases as the angular momentum increases beyond I = A + //. But this 
contradicts the fact that / = A + /i is a solution to the system when the function A has its 
maximum allowed value A max = 1. Thus, the maximum value of the angular momentum is 
X + fji. m 



Solution III. 2 When I = A + /x, there is a unique solution to ( \3. (\ )-( p7y[ ) with three possible 
cases for the deformation: (a) A > /i ; 7 = 7r/3 ; /3 = (A — /x)/3 and noncollective oblate 
rotation; (b) A < /i ; 7 = 0, /3 = (/i — A)/3 and noncollective prolate rotation; (c) A = /i, 
/3 = 0, 7 and tae rotation axis for the sphere are undetermined. 



From Eq. ( |3.10| ) for the quadratic Casimir, the deformation (3 and the total angular 
momentum I are evidently bounded. For no rotation, 1 = 0, the deformation attains its 
maximum /? max = + Xfj, + /j 2 /3. Although /3 = is a solution to the quadratic Casimir 

equation, it is not generally a solution to the cubic Casimir, Eq. ( p.lip . The exception is 
A = /i. Thus, the cubic Casimir invariant is essential to the upper bound for the angular 
momentum and the lower bound for the deformation. As proven above, the maximum 
angular momentum I = X + /i corresponds to the minimum deformation /3 m ; n = |A — /i|/3. 
In Fig. 0, p is plotted versus /; the solution curve is an ellipse. 



B. Analytical Solutions for Rotation Around One of the Principal Axes 

The system (|3 . Cj|) (|3T5|) can be solved analytically for the case of general (A, /x) and rota- 
tion around one of the principal axes, say the 1-axis. Assume that the angular momentum 
is directed along the principal 1-axis, I2 = I3 = 0, 1% = I. In this section, in contrast to the 
prior convention, the quadrupole moments of the principal axis solutions are not ordered. 
Thus the rotation axis is fixed, but q± may correspond to the short, long, or middle length 
axis. 

Solution III. 3 For rotations about the principal 1-axis, there are three analytical solutions 
of the algebraic system ( \3. 6j )-( \37Q ): 



10 



qi = \ — , q23 = — \ — ± -V A 2 — J 2 , < 7 < A, short axis rotation, (3.18) 

3 6 2 



2A + /i 2A + /i 1 i—z 

q 1 = H — , g 2 ,3 = 7, — ±-v/i -i , < J < /i, long axis rotation, (3.19) 

3 o 2 

A ~ A ~ M ± ^-y/(X + fi) 2 - P, < 7 < A + fi, middle axis rotation. (3.20) 



Ql = 7T~ , 92,3 



6 2 



Note that, since gi is the deformation along the rotation axis - short, long, or middle - the 
formulae for qi and q^ can be written as 



^±\VlL^P. (3.21) 



92,3 2 2 

It can be seen immediately that: 

1. The quadrupole moment along the rotation axis remains constant. 

2. At the maximal 7, viz. 7 max = A, 7 max = /i, 7 max = A + /i, these solutions are non- 
collective spheroidal, i.e. axially symmetric with respect to the axis of rotation. 

3. The three solutions correspond to bands rotating about the short, long, and middle 
axes of the 7 = nonrotating solution. For A > /i, the yrast band is described by 
the system rotating collectively about its short axis, < 7 < A. When the system 
rotates collectively about its long axis, the densities correspond to the bandheads, 
< 7 = K < /i. It is unclear if the densities describing rotation about the middle axis 
are found in nature; recall that such rotations are unstable in classical mechanics. 

4. Triaxiality is indicated by the 7 parameter. In Fig. the triaxiality is plotted versus 
the total angular momentum for short and long axes rotation in the (A, /j,) = (8, 4) 
irrep. Each band terminates in a phase transition to a noncollective rotation. For 
short and long axis rotations, the general expression for triaxiality is 



tan 7 



X + 2f,-V^P 

V 3 , short 

\ + 2ii + 3VX^P (3 22) 



- v y - 7 2 
\\ - , long. 
2A + ji 



5. When fi = 0, C 2 = f A 2 , C3 = |A 3 , < 7 < A, and the solutions scale as A: 




(3.23) 



When 7 = 0, the nucleus is a prolate spheroid. For small values of the angular 
momentum, / C A, the rotation is approximately collective prolate, but the shape 
is slightly triaxial. When the nucleus rotates with the maximum allowed angular 
momentum, 7 max = A, it is an oblate spheroid rotating around its symmetry axis (non- 
collective oblate rotation). The results for the deformation as a function of the angular 
momentum are summarized in Table |. 
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6. When A = 0, C 2 = f/i 2 , C 3 = — |/i 3 , < I < //, and the solutions scale as /i: 




(3.24) 



The initial shape (at / = 0) is an oblate spheroid. For small values of the angular 
momentum, I <C fi, the rotation is approximately collective oblate, but the shape 
is slightly triaxial. When / reaches its maximum value, J max = fi, the nucleus is a 
prolate spheroid rotating around its symmetry axis (non-collective prolate rotation). 
The deformation as a function of the angular momentum is given in Table [n[ 

7. When A = fi, C 2 = 2A 2 , C3 = 0. For a given angular momentum /, when (51, q 2 , qs) is 
a solution, so is (—qi, —q2, — 93)- 

• qi = and < / < 2A 

q 2 = iiv^A 2 - P 

q 3 = -q 2 (3.25) 



Here, the nucleus begins as a triaxial shape rotating around the middle axis. At 
/ max = 2A the nucleus turns into a sphere. 

• qi = A and < / < A. 



92 = __ + _VA^ 

Q^~-\V^P (3.26) 

In this case the nucleus is triaxial at / = and begins rotating around its long 
axis. At the maximum value of the angular momentum J max = A the nucleus is a 
prolate spheroid rotating non- collectively. 



C. Spheroidal solutions 



Analytical solutions can be found for spheroids either rotating about a principal axis or 
in a principal plane. Denote the two equal quadrupole moments by q and the unequal axis 
moment by —2q. The quadratic Casimir Eq. ( |3.7]) determines the deformation as a function 
of the angular momentum 

6q 2 + h 2 = C 2 . (3.27) 
Let K denote the component of the angular momentum along the symmetry axis. 
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1. Spheroids rotating around one principal axis 



Solution III. 4 K — 0: Collective rotation perpendicular to the symmetry axis. 



A solution must satisfy the cubic Casimir equation (3.8) 

- 24g 3 - 3ql 2 = 4C 3 . (3.28) 

Eliminating the angular momentum / from Eqs. ( |3.27D and ( |3.28| ) yields a single equation 
for the deformation, 6g 3 — 3C2<? ~ 2C 3 = 0. There is only one physical solution: 



I 2 = AXfi. (3.29) 

o 

When A > /i, this is a prolate spheroid; for A < fi, the solution is an oblate spheroid. The 
other two solutions to the system (|3.27| , |3.28|) , q = ~ (2A + /i) and q = — ~ (A + 2/i), are 
unphysical because they imply negative I 2 . 

Solution III. 5 K = I: Noncollective rotation about the symmetry axis. 

A solution satisfies the cubic Casimir equation ( |3.8| ) 

- 12g 3 + 3ql 2 = 2C 3 . (3.30) 

Eliminating the angular momentum / from Eqs. ( |3.27D and ( |3.30| ) yields a single equation 
for the deformation, 24g 3 — 3C2q + C3 = 0. All three solutions here are physical: 

g=(A + 2/i)/6, J = A; (3.31) 
q= -(2A + //)/6, I = fi; (3.32) 
q=(\- fJl )/6, I = \ + ^. (3.33) 

When q > the solution is an oblate spheroid; when q < 0, it is prolate. Note that these 
noncollective spheroids are the band terminations of the principal axis solutions ( [3.18D , 
flog ), and (ET2CD . 



Tilted rotation of spheroids 
Solution III. 6 zs between an<i /, z.e., t/ie rotation is tilted. 

Using the cubic Casimir equation, Eq. (|3.8|), the squared projection K 2 may be solved for 
function of the angular momentum /, 

= P + 4Cj 

9 9q K ' 



When A > /i, a prolate spheroid (q < 0) solution exists in the interval pi < I < v^A/x with 
the squared projection of the angular momentum on the symmetry axis 
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K 2 = /x 2 + -(I 2 - /i 2 ) + ^(A - /i)(A + 2/x) 



1 - 



2A + /i 



v /(2A + / i)2-3(J 2 -/i 2 ) 



/i 2 + 



(2A + /i) 2 



(A-/x)(A + 2/x) ,_ 2 



(2A + n)< 



.2\2 



(3.35) 
(3.36) 



For I = fi, the prolate spheroid rotates noncollectively about its symmetry axis, Eq. fl3.32| ) . 
At a maximum / = \/4A/x, the prolate spheroid is rotating perpendicularly to its symmetry 
axis. In Fig. [5] the projection if is plotted versus I for the (A, /x) = (8, 4) prolate spheroids; in 
Table [TTT| the projection and the deformation are given for this case. Note that the projection 
is approximately constant K pa 4 until I nears the top of the band; this is evident from the 
Taylor expansion (|3.36|) . Thus there is a K ~ /x band of prolate states for /x < 7 < 2-y/A/x. 
The band may be viewed from the top down: Start with a collective rotation of a prolate 
spheroid with I = 2-y/A/x- As the body rotates more slowly, it rapidly acquires a component 
K pa 4 along the symmetry axis. As the magnitude of the angular momentum becomes 
smaller, its direction becomes more aligned with the symmetry axis. The band terminates 
in a noncollective rotation when I = K = /x. 

When A > /x there are also tilted oblate solutions for A < J < A + /x. The two endpoints 
correspond to noncollective rotation and termination of the short and middle axes bands 
O and (|3~33D . 

When A < /x, an oblate spheroid (q > 0) rotates with the projection 



A 2 + -(/ 2 -A 2 )-^-(A-/x)(/x + 2A) 
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A 2 + 



A 2 



27 
•(/ 2 - V 



2/x + A 



, (A //)(// • 2 A) /; . x2i2 
(2/x + A) 2 ^ ' W (2/x + A) 4 



^ (2/x + A) 2 - 3(/ 2 - A 2 ) 
(J 2 - A 2 ) ; 



(3.37) 
(3.38) 



For I = X, the solution is an oblate spheroid rotating noncollectively about its symmetry 
axis. At a maximum / = \/4A/x, the oblate spheroid is rotating perpendicularly to its 
symmetry axis. When A < / \/4A/x, is approximately A. From the Taylor expansion 
( |3.38| ), the projection is approximately constant K pa A until one nears the top of the band. 
If A = fi, the cubic Casimir is zero, and 



K 2 



I 2 - u 2 

/x + 
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(3.39) 



As the angular momentum varies from a minimum I — /x to a maximum I = 2/x, the 
projection of the angular momentum on the symmetry axis varies from K = /x to K = |/x. 



IV. ENERGY FUNCTIONAL I 

The particular solutions enumerated in § [TIT] correspond to simple kinematical situations. 
The actual physical densities must be determined from a dynamical argument. An equilib- 
rium density for a rotating body in the su(3) model is a critical point of the energy functional 
E[p] on the surface of admissible densities. For a body with constant moments of inertia, 
the energy is 
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E[p] = AJl + A 2 I 2 + A 3 I 2 3 , 



(4.1) 



where A±, A2, A3 are real constants. Suppose the total angular momentum is J, the admis- 
sible densities lie on the surface determined by the quadratic and cubic Casimirs Ck, and p 
is a critical point of E[p] on the surface. By the Lagrange multiplier theorem, there are four 
real constants a, b, c, d associated with each of the constraint equations (|3.6|-|3~^) such that 
p is a critical point of the functional 



F[p] = E[p] +a( qi + q 2 + q 3 ) + blj2 

V k 



+d (I 2 + II + If - I 2 ) (4.2) 

with a free variation on (gi, q 2 , q 3 , I\, I 2 , h)- Altogether there are 10 unknowns, including 
the Lagrange multipliers, that must satisfy the system of 10 equations, 

OF f 3 \ 

° = dq~k = a + 2bqk + ° \ q2k " 4 J V (4 ' 3) 

= f£ = 2/ fc (A k - -cq k + d) (4.4) 



dl k y 4 

plus the four constraint equations. Although the number of equations in the system is large, 
each is just a polynomial of low degree in the variables. Thus analytic solutions may be 
determined rather easily in many cases. 



Indeed each of the kinematical solutions enumerated in § [TTT] are particular solutions to 
the Lagrange multiplier system: 



1. Rotations about a principal axis, say the 1-axis, given by Eqs. (|3.18| , |3.19| , |3.20|) are 
critical points of (|4.2|) when d = Ai and a = b = c = 0. 



2. Spheroidal solutions, g 2 = q 3 = q 7^ and q\ = —2q, are obtained when a = b = c = 
and 

h{Ai ~d)= h(A 2 -d) = I 3 (A 3 - d) = 0. (4.5) 

Principal axis rotation is one type of solution. When two energy parameters are equal, 
say A\ = A3, tilted rotations are also critical points for d = A\ and I 2 = 0. 

It is interesting to compare the classical rigid rotor to an su(3) rotor. A classical rigid ro- 
tor in equilibrium is a critical point of the energy ( |4.1|) subject to the constraint 7| = I 2 . 
This is mathematically equivalent to the su(3) constraint problem when a = b = c = 0, 
viz., Eq. ( |4.5|) . The su(3) problem allows for additional critical points. The mathematical 
origin of the extra solutions is the different constraint system for the admissible su(3) sur- 
faces. The physical origin is related ultimately to the shell model, which is responsible for 
the su(3) quadrupole operator (|2.1|) instead of the major shell-mixing quadrupole operator 
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Suppose the kinematical solutions corresponding toa = b = c = are excluded from 
consideration. In this solution to Eq. ([O) requires 



= 4 (q 2 - qi ) (q 3 - qi ) (q 3 - q 2 ) (4.6) 
+A 2 (92 - Qs) + if (q 3 - qi) + if (?i - <? 2 ). 



In conjunction with Eqs. (|3.8| , p.9|) this condition enables the determination of the body-fixed 
angular momentum components in terms of the quadrupole moments: 

T 2 = P , 4 (gf + 2q 2 ql + 2 g| g 3 + g 2 4 ) - 2 gi C 3 
1 3 + 3 (ql + q 2 q 3 + q 2 ) 

and similarly for I 2 ,I 3 by cyclic permutation of 1,2,3. After eliminating the Lagrange 
multipliers the system reduces to a single independent equation, 

= I 2 {{A 2 - A 3 ) qi + (A 3 - A 1 )q 2 + {A 1 - A 2 )q 3 ) , (4.8) 

where 7i, /3 are assumed nonzero. There are two cases to consider as one of the two factors 
in the above equation must vanish. 

Solution IV. 1 Principal plane rotations of triaxial bodies. One of the components of the 
angular momentum vanishes, say the 2— axis projection, I 2 = 0. Every critical point of the 
principal plane system is determined by a real root q% of the eighth degree polynomial 

= 576 q^ + 288 (i 2 - 2 C 2 ) q\ + 192 C 3 q\ + (36 1 4 - 168 C 2 1 2 + 192 C 2 2 ) q\ + (72 C 3 I 2 - 144 C 2 C 3 
+ (-12C 2 I 4 + A8C 2 I 2 + 16 C 2 3 - 48 Cf) q\ + (8 C 2 C 3 -AC 2 C 3 1 2 ) q x + C 2 1 4 + (4 C\ - 4 C 2 3 ) / 
-8C 2 C 2 3 +AC 4 . ( 

When q% is a real root, the other two quadrupole moments and angular momentum compo- 
nents are given by 

(I 2 -2C 2 ) (Qq 2 -C 2 ) 

q2 = 4(6g3 + C 3 ) ' * = ^ 10 ) 

j2 = ( q2 - 9 3 ) I 2 - 4 [ql ( q 2 - gg) + <j{q z - gi) + qj ( <h - g 2 ) ] j2 = p _ p u u \ 

3 3^2 

A valid solution is attained whenever the squares of the angular momentum components 
If, I 2 are nonnegative. 

Since I 2 = there is a reflection symmetry: when q\ is a valid solution, q 3 also satisfies 
the eighth degree polynomial. Thus, solutions come in pairs except in the spheroidal case 
qx = q 3 . To eliminate this trivial redundancy, we choose q% < q 3 . In Table |TV| two sequences 
of solutions are given for the case (A,/i) = (8,4). There are no valid solutions for angular 
momentum I less than fi = 4. When 1 = 4 the noncollective prolate solution ( p.32| ) begins 
a band of principal plane triaxial solutions that terminates at / ~ 10.08. When I = 8 a 
second band emerges starting from the noncollective oblate state ( |3.31| ) and also terminating 
at / ~ 10.08. Above the critical angular momentum / w 10.08, the pair of real roots turns 
into a pair of complex conjugate roots. 
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There is a third band of principal plane triaxial rotors that begins at / = y/AXfi which is 
the termination of the sequence of prolate spheroidal solutions ( p. 29 ). This band terminates 



at / = A + /i in the noncollective oblate state. In Figure f| the projection I\ of the angular 
momentum on the short axis for this band is plotted versus the angular momentum I for the 
case (A, fi) = (8, 4). In Figure |3] this band is shown as the continuation of the prolate spheroid 
band to triaxial principal plane rotors when / > y/AXJI; I 2 — K — is the projection of the 
angular momentum on the long axis. 

Solution IV. 2 Constant-^ bands. If all three components of the angular momentum are 
nonzero, Ik ^ 0, then 

(A 2 - A 3 ) gi + (A 3 - A l )q 2 + {A 1 - A 2 )q 3 = 0. (4.12) 

When two of the energy coefficients are equal, say A 2 = A 3 , this equation forces a spheroidal 
solution, q 2 — q 3 . When the coefficients are all distinct, the system may be solved for the 
quadrupole moments: Fori,j,k cyclic, 

q i = (2A i -A j -A k )u, (4.13) 



u 2 = 2C 2 -I 2 

12 (Al + A% + Al-A 1 A 2 -A 2 A 3 -A 3 A 1 y { ' ] 

Since the ratios among the quadrupole moments are constant, the triaxiality parameter 
7 is fixed and independent of the angular momentum. The angular momentum components 
in the intrinsic frame are determined from Eq. ( |4.7| ); a valid solution demands that the 
square of each angular momentum component is nonnegative. These fixed-7 bands typically 
start at rather large angular momentum and terminate at the maximum angular momentum 
I = X + fi. 



V. ENERGY FUNCTIONAL II 



A second natural choice for the energy functional is a polynomial in the so (3) integrity 
basis of su(3). By definition, any rotationally invariant polynomial in the su(3) enveloping 
algebra is a function of the so (3) integrity basis. The basis consists of five functionally- 
independent polynomials: the two Casimirs, C 2 ,C 3 , the square of the total angular momen- 
tum L 2 , and two scalars of degrees three and four [?], 

X 3 = L i Q^>L j (5.1) 
X 4 = L,Qg ) QgL fc . (5.2) 

The eigenvalues of the hermitian operators A3 and A4 in a representation are not integers 
[?]. On a level surface the Casimirs are constants and the other three integrity basis members 
may be evaluated in the principal axis frame, L 2 = I 2 , X 3 = q k I\, A 4 = q\l\. A simple 
energy functional used in nuclear structure applications is a linear combination of J 2 , A 3 , 
and A 4 , 
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E[p] = AI 2 + BX 3 + CX 4 , 



(5.3) 



where A, B,C are real constants [?]. Note that the moments of inertia for this energy 
functional are not constant, and they depend on the deformation. When C ^ 0, the energy 
may be expressed conveniently as 

E\p\ = AI 2 + C{uX 3 + X 4 ), (5.4) 

where v — B/C. Since I 2 is constant on the constraint surface and C determines an energy 
scaling, the critical points depend only on the real parameter v. The critical points of Q5.4| ) 
on the constraint hypersurface (|3.6| - 13.9| ) in the dual space are determined by the Lagrange 
multiplier theorem in a way similar to Eq. (|4.2p . 

Once again, the kinematical solutions of § |TTTj are critical points of ( |5.4j ). 



1. Rotations about a principal axis, viz., Eqs. ( |3.18| , |3.19| , p.20|) , are critical points for the 



multipliers, a = 3cq 2 q 3 , b = 3cgi/2, d = q\(q\ + v + 3c/4), and 

I 2 (2 qi + v) 
C 3(2q 2 + q 2 q 3 -P/AY 



(5.5) 



When the denominator in the above expression for c is zero, there is still a principal 
axis solution, but only when the parameter v = —1q\. For example, the band given 
by Eq. (gTg) for < / < A yields 



9/i(A + /i) 

When [i = the energy parameter must be v = 2 A/3. 

2. Spheroidal solutions are critical points only when I 2 = I3. In this case the Lagrange 
multipliers may be chosen to yield the spheroidal solutions. 

The equations for general triaxial rotors are obtained after eliminating the Lagrange 
multipliers from the system. Ignoring the special case of spheroidal solutions, a solution 
requires that the rotation is in a principal plane. 

Solution V.l Principal plane rotations. When I2 = the body-fixed angular momentum 
components 1%, I3 are determined from Eqs. In addition to the constraint equations 

({Of, \3. the quadrupole moments satisfy 



= 12(qf + ql - 4 qi q 2 q 3 ) - 3I 2 q 2 

+Qv(q 3 -q 2 )(q 1 -q 2 ) + 2C 3 . (5.7) 

Critical points of the energy functional are shown in Figures fl, ||, [7| for (A, fi) = (8, 4) and 
for three values of v. 20/3, 0, and —4/3. For each chosen v there are two plots in the figure: 
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one of the triaxiality parameter 7 versus the angular momentum, and the other of the square 
of the ratio versus the angular momentum. 

In these figures only the principal plane Solutions V.l are drawn. There are, of course, 
the usual kinematical bands for principal axis rotations and spheroids which are not shown. 
The equations are symmetrical under a simultaneous interchange of q l and q 3 , I x and I 3 . 
As a result, the plots of (Ii/I) 2 are reflection-symmetric with respect to the (Ii/I) 2 = 0.5 
horizontal line. One of the curves may be regarded as a plot for (Ii/I) 2 while its reflection 
is a plot of (/3/J) 2 . When (Ii/I) 2 = 0.5 the tilting of the rotation axis in the 1 — 3 plane is 
a maximum; when — or 1, the rotation is about a principal axis. 

If v = 20/3, there are two bands, Figure |5|. The first band begins at the triaxial 1 = 
state and ends at / = 8 as a noncollective oblate spheroid. Note that maximal tilting 
is attained for the 1 = density - in vivid contrast to the principal axis solution. The 
second band starts at I = 4 as a noncollective prolate spheroid and ends at I = 12 as a 
noncollective oblate spheroid. Although the two end points of this band are principal axis 
rotations, the densities in between are tilted as the direction of the angular momentum 
vector varies continuously from alignment with the long axis (J = 4) to alignment with the 
short axis (J = 12). 

When v = —4/3 there are two bands, but there is a gap between them. The first band 
begins as a triaxial maximally-tilted rotor at / = and ends as a noncollective prolate 
spheroid at I = 4. There is a sharp phase transition to the I = 4 principal axis rotor from 
a highly tilted configuration. The second band is oblate-like from I = 8 to / = 12 while the 
rotation is slightly tilted. 

If v = there are two bands, no solutions below 1 = 4, and a gap between the bands. 
The first band begins as a triaxial, maximally-tilted rotor at I = 4, increases in angular 
momentum to I ~ 6 before falling back to the principal axis noncollective prolate solution 
at I = 4. The second band begins at / = 8 and ends at / = 12; these oblate-like state are 
only partially tilted. 



VI. COMPARISON WITH REPRESENTATION THEORY 

The correspondence between density matrix theory and the irreducible representations of 
su(3) is quite close, but not perfect. In this section the similarities and differences between 
the two are discussed. 

Among the continuous family of surfaces of admissible densities, a representation may 
be associated only with those that satisfy a generalized Bohr-Sommerfeld quantization rule 



[?]. Consider the diagonal density matrix P(A At ), Eq. Q2.12|) , which is contained uniquely in 
each surface. For each group element in the subgroup of diagonal matrices 

g = dmg(e ie \e ie \e ie *), det<? = l, (6.1) 

where #1, 82, #3 are real numbers, define the complex number 

X M {g) = exp {i(2\ + ii)9x/3 + + //)0 2 /3 

-i(X + 2/^/3} . (6.2) 

The quantization condition is that \^ (f?) is a character of the diagonal subgroup, i.e., 
whenever g is the identity element of the group, x (<?) must equal one. When g = 
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diag(e j6> , e~ ie , 1) is a closed circle, < 6 < 2n, the character x (9) = e%Xd must equal 
one at 9 = 2tt, or A is an integer. Similar reasoning for g = diag(l, e ld , e~ ld ) shows that fi 
must be an integer. A level surface of the Casimir invariants satisfies the Bohr-Sommerfeld 
quantization rules if and only if A and fi are nonnegative integers. Hence the surfaces of 
admissible densities that satisfy the quantization condition are in one-to-one correspondence 
with the irreps. 

The irreducible representations of the compact algebra su(3) are finite dimensional; 
hence, the spectrum of any hermitian operator acting on the representation space is bounded. 
On a surface of admissible densities, the range of any real-valued continuous function is 
bounded from above and below, because these surfaces are closed and bounded manifolds. 
Thus the values of physical 5u(3) observables have a bounded range, whether hermitian 
operators acting on irreducible representation spaces or real-valued functions on surfaces of 
admissible densities. In particular, the range of the total angular momentum is the same 
bounded interval, < / < A + /i, in an irrep of su(3) and on the corresponding level surface 
of the Casimirs in the dual space su(3)*. 

The differences between representation theory and density matrix theory are more ev- 
ident from a comparative analysis of the angular momentum decompositions. For both 
su(3) irreducible representations and the adiabatic rotational model, the decomposition is 
achieved by angular momentum projection from a fixed intrinsic state. When the intrinsic 
state has a sharp quadrupole deformation, as it does in the adiabatic rotational model, the 
projected states \IKM) are orthogonal. In addition, each projected state is an eigenstate, 
belonging to the eigenvalue K 2 , of the hermitian operator I 2 that corresponds to the square 
of the third component of the body- fixed angular momentum [?], 

(<?3 - <?i)(<?3 - 9 2 )4 2 = X 4 + q 3 X 3 + q x q 2 I 2 , (6.3) 

where the quadrupole moments (q x , q 2 ,(lz) are rea ^ constants determined by the intrinsic 
state. 

In contrast, angular momentum projection from a su(3) highest weight vector yields a 
nonorthogonal basis that requires Gram-Schmidt orthogonalization [?] . Moreover K 2 is not 
an eigenvalue of J| which, indeed, is no longer well-defined for su(3). These difficulties arise 
because the 5u(3) intrinsic state, the highest weight vector in an irrep, does not have sharp 
values for the quadrupole moments. However, when a vector's angular momentum / is small 
compared to A and/or /i, its deformation is little changed from that of the highest weight 
vector, Eq. (|3.7|). Thus, for / <C max(A,/i), the 5u(3) projected states are approximately 
orthogonal and K 2 is an approximately good quantum number. According to Elliott [?], 
the possible values of K in an su(3) irrep (A, /i) are 

K = min(A, /i), min(A, ji) — 2, 1 or 0. (6.4) 

The sequence of angular momentum states associated with a K-band is: 

K = : I = max(A, //), max(A, //) — 2, 1 or , fi . 

K^O : I = K,K + 1,...,K + max(A,/i). ^ ' 

The density matrix theory shares some properties with both the geometrical rotational 
model and su(3) irreducible representations. The densities have well-defined values for the 
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deformation, but, on a level surface of the su(3) Casimirs, the deformation depends on the 
angular momentum. In su(3) density matrix theory, Eq. Q6.31 ) for J| is an identity when the 
operators are replaced by their corresponding functions on the dual space and, instead of 
remaining constant, the quadrupole moments now vary. For small angular momentum, as 
long as the deformation is constant, the geometrical model, the density matrix theory, and 
the irreducible representations of su(3) are equivalent theoretical descriptions. For larger 
angular momentum, the connection between the 5u(3) and geometrical models breaks down. 
Yet, because the deformation changes on a surface of admissible densities, the matrix theory 
maintains its correspondence with su(3) irreps. 

Suppose A > In the density matrix theory, the band heads < I — K < jj, are densities 
rotating about the long principal axis, Eq. ( |3.19|) ; the last bandhead is at I = K = // when 



the body is a prolate spheroid rotating about its symmetry axis. The band heads correspond 
to the Elliott labeling of K bands, Eq. fl6.4T) 

The Elliott K — band corresponds to the density matrices which describe bodies 
rotating about the short principal axis, Eq. ( p.l8|) ; this band terminates at I — A when an 
oblate spheroid rotates about its symmetry axis, Eq. ( p.5| ). 

For K ^ bands, the relationship between density matrix and representation theory is 
not as clear. An obvious candidate for the Elliott K — \i band is the sequence of prolate 
spheroid densities which attains a maximum / = y/A\jx, § |III C 2\ At this maximal / 



the body rotates collectively about its short axis. Because the projection of the angular 
momentum on the symmetry axis does not maintain a constant value K = fi, Figure |3|, 
the geometrical interpretation as a fixed K band is not supported. In addition the density 
matrix band terminates before I = X + fi. In § [TV] a sequence of triaxial principal plane 
rotational densities which starts at / = -\/4A/i and terminates at / = A + /i was found. In 
our view the density matrix theory indicates that the Elliott K — \i band in a su(3) irrep 
does not exist in the strict geometrical model sense, although it is a useful concept when 
H < I < V4A/2. 

For the other K bands, there are more discrepancies. In particular, K bands in the 
density matrix theory need not terminate at / = K + A. However, when A is large compared 
to both /j, and the angular momentum /, approximate analytic solutions to the fundamental 
system, Eqs. (|3.6| - |3.9| ), can be found that describe tilted rotation of triaxial rotors. For 
X ^> I > K, there are approximate solutions for which the quadrupole moments in the 
principal axis frame are 
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2 A + n I 2 — K 2 



3 2(2A + /i) 



ygHZ! + (6 . 6) 

H 6 2 4 (2 A + n) K ' 

2 A + a Jfi 2 - K 2 I 2 - K 2 
<?3 = 7, ^ h 



6 2 4 (2A + /x) 

and the body-fixed projections of the angular momentum are 

1 (2X + rf K ' 
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l\ = \[l 2 - K 2 ) (l ) (6.7) 

2 V ' \ (2A + /i) 2 ' 



p 



H 2 



(2X + nY 



The quadratic Casimir equation is satisfied to o(A -2 ) and the cubic Casimir equation to 
o(A -3 ). For A large compared to / and /i, these triaxial rotors are approximately prolate 
spheroids. The solutions form a band in which the component of the angular momentum 
along the long near-symmetry axis is approximately K. 

To distinguish among basis vectors with the same total angular momentum in an su(3) 
irrep, states may be chosen to be simultaneously eigenstates of some element of the su(3) D 
so (3) integrity basis, e.g., X 3 , X 4 or some more complicated function of the integrity basis [?]. 
Diagonalizing just X 3 is the simplest choice [?]. For axially symmetric states the hermitian 
operator Jf is related to X 3 [?], 

I 2 = 1 -( 1 - det QV)- -V3 X 3 + h 2 . (6.8) 

A stringent test for the density matrix method is to compare the eigenvalues of X% 
calculated in representation theory with the simple geometrical densities. In Table |Vj, the 
smallest eigenvalues of X 3 for each angular momentum state from / = to / = 8 in the 
(8,4) irrep are compared to the X 3 values for K = geometrical densities (|3.18|). Since X 3 



is a cubic polynomial in the enveloping algebra and not an element of the algebra itself, the 
evaluation <p(X 3 ) of this function in the dual space must be defined now. These values are 
calculated for rotation about the short principal 1-axis by 

<p(X 3 ) = qj 2 . (6.9) 

The error in the density matrix calculation rises to a maximum of just 10% for the 1 = 8 
state. 

For the I = 12 state of the (8, 4) irrep, the eigenvalue of X 3 is —234; for the noncollective 
rotation of the oblate spheroid density, (p(X s ) = (/i — A) (A + /i) 2 /3 = —192. 

For the bandheads J = K = 0, 2, 4, the states in the representation space are taken 
to be the eigenvectors belonging to the maximal eigenvalues of X 3 , 0, 25, 108, respectively. 
For the corresponding geometrical densities given by rotation about the long principal axis, 
(p(X 3 ) = 0, 27, 107, as calculated via the analytic formula ( 3.1 9| ) . The error is negligible for 
small angular momentum. 

In the representation space the U K = 4" band is chosen to be the sequence of angular 
momentum states belonging to the maximal eigenvalue of X% for each angular momentum. 
In Table |VT|, these eigenvalues are compared to the values of (p(X 3 ) for the sequence of tilted 
prolate spheroids, 

V(X 3 ) = q (I 2 - 3K 2 ) , (6.10) 

where q is the moment for the equal short axes and K is the component of the angular 
momentum along the long symmetry axis. 

The U K = 2" band in the representation space is taken as the sequence of states with 



intermediate values for the X 3 eigenvalue, see Table VII . In this case the geometrical 
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densities are assumed to be principal plane triaxial rotors, I2 = and I3 — K — 2 where 
the 3-axis is the long axis. The complex values at the end of the K = 2 band for the values 
of the quadrupole matrix are due to the restriction of the rotation in one plane. They show 
that the assumption I 2 = is not correct. There are non-zero projections of the angular 
momentum on all three intrinsic axes. 

When fi = 0, e.g., 20 Ne for which (8,0) is the dominant irrep, the angular momentum 
is multiplicity-free, I = 0, 2, . . . , A. In this case, Bargmann and Moshinsky [?,?] have given 
an analytic formula for the X 3 eigenvalues: X 3 = (A/3 + 1/2)1(1 + 1). The density matrix 
approximation is X 3 = qil 2 = (A/3)/ 2 . The difference is due to the omission of commutator 
terms of lower degree in the density matrix approximation - I 2 instead of the quantum 
/(/ + 1), and instead of 1/2. 

VII. COMPARISON WITH CRANKED ANISOTROPIC OSCILLATOR 

In conventional mean field theory nuclear rotational motion is modeled in a simple way 
by cranking the anisotropic harmonic oscillator around one axis. For rotation with constant 
angular velocity uj around the 1-axis, the Hamiltonian (Routhian) for one nucleon with mass 
m in the rotating frame is 



where u>k (k — 1, 2, 3) are the oscillator frequencies. For many fermions the model wavefunc- 
tion is a Slater determinant given by occupying the orbitals of h u , and it is an eigenstate 
of the one-body operator H w = ^ a h u , where the sum is over the particles. When the de- 
formation is not too large, mixing between major oscillator shells can be ignored [?,?], and, 
as a function of the total angular momentum /, the energy of a system of many nucleons in 
the anisotropic potential simplifies to 



where denotes the sum of the quanta (rifc+1/2) over all occupied orbitals. J max = |S 2 — S 3 | 
is the maximum angular momentum of the rotational band, < / < J max . Applying 
Feynman's lemma, 



the expectations of the dimensionless quadrupole moment in the rotating frame are 






(7.2) 




(7.3) 




(7.4) 



(7.5) 



(7.6) 
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These expectations are exactly the values of the quadrupole moment found with the density 
matrix method, Eqs. ( p. 18 , p. 19 , |3.20|) . For example, when Si > X 2 > S 3 , the rotation 



is about the long axis, Eq. ( |3.19| ), where A = £1 — E 2 , n = £ 2 — S3 = I max - But their 



derivations involve different assumptions. The su(3) density derivation shows that principal 
axis solutions are an immediate consequence of kinematics and the restriction to su(3) 
admissible densities. The constraints imposed by the quadratic and the cubic Casimirs are 
essential. The principal axis solutions are critical points for both energy functionals I and 
II. The natural inference is that the principal axis solutions should be critical points for 
any physically reasonable energy functional. In the density matrix theory the energy is 
proportional to I 2 for both functional I and II, 

Ej(I) = qj\ Ejj(I) = (A + B qi + Cq 2 ) I 2 . (7.7) 

In contrast the anisotropic oscillator derivation is based on a specific assumption about 
the energy. When major shell mixing is ignored, the cranked anisotropic oscillator is com- 
patible with su(3) dynamical symmetry. By making a canonical (symplectic group) trans- 
formation and expressing the su(3) generators in the transformed basis, the Hamiltonian 
of the cranked anisotropic oscillator (CAO) becomes an element of the algebra. Hence the 
su(3) density prediction for the quadrupole deformation is expected. Note, though, that the 
energy of the cranked anisotropic oscillator is not proportional to I 2 . When the oscillator 
frequencies are optimized, the self-consistent energy of the cranked anisotropic oscillator is 



1 



1/3 



E CAO = 3Hw (Si (£ 2 £ 3 + j , (7.8) 

where uJq = ujiuj 2 ujz is constant. 

When the deformation is large the assumption of no major-shell mixing is not valid. The 
anisotropic oscillator energy and the expectations of the angular momentum and quadrupole 
moment can be evaluated analytically in terms of the angular velocity uo [?]. The final 
results cannot be expressed as analytic functions of the angular momentum, but numerical 
calculations are elementary. For rotation around the 1-axis, the numerical values for the 
quadrupole deformations from density matrix theory, Eqs. ( p,18|J3.1^|J3.2(J| ), and from the 
exact problem of the cranked anisotropic harmonic oscillator are compared for two even- 
even cis-shell nuclei: 20 Ne and 24 Mg. 

The case of the highly deformed nucleus 20 Ne is presented in Table |VTII]. This nucleus 



is described by the (8, 0) su(3) representation, which, for collective rotation around the 
short 1-axis, corresponds to £1 = S 2 = 14, £ 3 = 22. The deformations from the density 
matrix method are calculated from Eq. fl3.18p . These are compared to the values calculated 



numerically with exact cranking of the anisotropic oscillator [?]. The agreement is excellent. 
In particular, the band end-points - prolate spheroid at I = and noncollective oblate 
spheroid at / = 8 - are in perfect agreement. 

The case of the tri axially deformed nucleus 24 Mg is more interesting. The dominant 
su(3) representation for this nucleus is (8,4), which corresponds to £ 3 = 28, S 2 = 20, 
Si = 16. In Table |IX| the deformations from the cranked anisotropic oscillator are compared 
to the values from Eq. ( p . 1 8[ ) for collective rotation around its short 1-axis. The differences 



between the two theories are negligible. 



24 



In Table [X] the deformations from the cranked anisotropic oscillator are compared to 
the values from Eq. ( p.!9| ). The agreement is excellent and shows that Eq. Q3.19D can be 
interpreted as the formula for the band heads: The nucleus rotates around its long axis. 
At / = 4 the nucleus is a prolate spheroid rotating non-collectively. This is the maximum 
angular momentum at which a band occurs for the dominant representation (8,4) according 
to the Elliott model. 

VIII. CONCLUSION 

The density matrix method provides a simple geometrical interpretation for the rota- 
tional states in 5u(3) irreducible representations. Each density in the dual space has a direct 
physical interpretation as the expectation of observables in the algebra. But the quantum 
superposition principle and state fluctuations are not incorporated directly into the density 
theory because the admissible densities do not form a vector space. Nevertheless, the den- 
sity formulation reproduces many properties of the quantized irreducible representations. 
This situation is similar to Hartree-Fock and its relationship to the quantized shell model. 
Although Hartree-Fock was founded on the independent fermion assumption, density matrix 
theory shows that, when viewed from an appropriate perspective, the essential character of 
the mean field method does not demand this assumption. In fact density matrix theory can 
be applied effectively to describe collective rotational states. 

SU (3) density theory is more tractable and than either the irreducible representation 
theory or the cranked anisotropic oscillator. The densities are solutions to a system of al- 
gebraic equations that are given immediately by the model's ansatz. In contrast, even for 
principal axis rotations, the cranked anisotropic oscillator requires an extended argument 
involving energy minimization to attain the same conclusions. The irreducible representa- 
tions, determined from the theory of highest weights, are difficult to work with in the non- 
canonical angular momentum basis. For example, the eigenvalues of X3 in an irreducible 
representation are difficult to compute, yet the values from density matrix theory are a 
back-of-the-envelope calculation. This illustrates the power of the density matrix theory 
and indicates its potential for analyzing more complex algebraic models. 

Density matrix theory transforms the angular momentum decomposition problem into a 
geometrical analysis of the range of the angular momentum function on admissible surfaces. 
In many cases, such as SU(3), this is technically easier than the mathematical procedure 
of identifying irreducible subspaces of 5*0(3) within an irreducible representation of the 
model's algebra g. The latter is a difficult task when 5*0(3) is not canonically embedded in 
0- 

In density matrix theory, the cubic su(3) Casimir plays a role equal to the quadratic 
Casimir in determining the set of admissible densities. In the usual su(3) model, the 
quadratic Casimir has a distinguished part because the quadrupole-quadrupole interaction 
is a linear combination of C2 and I 2 . 

The density theory allows for many solutions, including tilted rotation in a principal 
plane and more intricate rotational configurations. The physical interpretation of these 
solutions is simple and may be adopted for the corresponding state vectors in irreducible 
representations. 
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The density matrix method may be applied to any dynamical symmetry algebra g. The 
admissible densities of the model are a level surface of the Casimirs in the algebra's dual 
space. For a semisimple Lie algebra, the dimension of the generic level surface O p equals 
the dimension of the algebra minus the rank of the algebra. There are also singular level 
surfaces whose dimension is even smaller, e.g., A = or /i = level surfaces are just four 
dimensional. In nuclear applications the angular momentum algebra so (3) is a subalgebra 
of q. Hence, a rotationally invariant energy functional enables a further reduction in the 
dimension to dim O p — 4 after rotation of the system to an intrinsic frame (a reduction by 
three, the dimension of the rotation group, and restriction to fixed angular momentum I). 
For the su(3) problem, the dimension of the algebra is eight, its rank is two, and the effective 
dimension is two for a rotational scalar energy functional. 

In future work the method will be applied to other algebras relevant to nuclear structure 
science, e.g. the general collective motion algebra gcm(3) corresponding to the extended 
Bohr-Mottelson model, which includes quantum vorticity [?], the symplectic algebra sp(3, R) 
[?,?], and the interacting boson model u(6) and its subalgebras u(5), 50 (6), and 50 (5) [?]. 
It should be emphasized that the only restriction is that the physically relevant degrees of 
freedom span a Lie algebra of observables. 

The Hohenberg-Kohn theorem of density functional theory was generalized to establish 
the existence of an energy functional for arbitrary dynamical algebras whose minimum is 
the density of the exact ground state [?]. But, like the original Hohenberg-Kohn result 
[?], this is an existence theorem for which an explicit construction of the density functional 
from the Hamiltonian is not known. Nevertheless, like the Hohenberg-Kohn theorem, it 
suggests a promising avenue of research to solve complex many-body problems. Two energy 
functionals for su(3) were considered in this paper, one motivated by the classical theory of 
rigid rotations, the other - from the mathematics of integrity basis theory. 

The surface of admissible states is, in fact, an orbit of the coadjoint action of the Lie 
group in the dual space [?,?]. These orbit surfaces are equipped with a Poisson bracket 
and symplectic structure. In particular, each surface is always even-dimensional and admits 
canonical coordinates. An energy functional defines a Hamiltonian function, and, from the 
Poisson bracket, a Hamiltonian dynamical system. Hence the dynamics of density matrices 
is well-defined. Even though this paper studies equilibrium densities in the su(3) theory, nor- 
mal mode and other dynamical properties may be investigated too. These are the analogues 
of the random phase approximation and time- dependent Hartree-Fock from conventional 
mean field theory. 

Finally the surfaces that satisfy the generalized Bohr-Sommerfeld quantization condition 
may be used to construct explicit irreducible representations. The procedure to obtain the 
irreps is called geometric quantization [?,?,?]. The method was applied in prior work to 
determine the irreducible representations of the rotational and Bohr-Mottelson theories [?]. 
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APPENDIX A: DUAL SPACE FUNCTIONS 



One important concept in density matrix theory is that of a coordinate function. For 
each Lie algebra element Z e su(3) there is a real- valued "coordinate" function (p(Z) defined 
on the dual space: the value of the function (p(Z) at the point p G su(3)* is defined by 

<p(Z)(p) = {p,Z) =ti (j>Z). (Al) 

These functions separate points, i.e. if p\ and p 2 are distinct densities, then there exists a 
Lie algebra element Z such that ip(Z)(pi) ^ ip(Z)(p 2 ). The physical interpretation of <p{Z) 
is that it is the real- valued function corresponding to the hermitian operator o{Z). The 
value of the observable Z at the density p is the number ip(Z)(p). 

Polynomials in the algebra generators are relevant for many physical applications, e.g., 
the Casimir invariants and the cubic and quartic scalars X 3 and X 4 . In this appendix the 
extension of ip from the domain of Lie algebra elements to the domain of polynomials is 
defined. 

First, the domain of ip is extended to the symmetric algebra S(g) of the Lie algebra g. 
S(g) is a commutative associative algebra with elements that are symmetric polynomials of 
elements of g (the order of multiplying the elements of g is not important) [?]. The extension 

ip : S{g) — >C~( fi r,R) (A2) 

is defined in the following way: if {Z ± , Z 2 , . . . , Z n } is a basis for the Lie algebra g and 
(ei, e 2 , • . . e n ) is a set of integers, then 

<p(Z?Z? . . . Z?) = p>{Z x )^p>{Z 2 T . . . <p(Z n y™. (A3) 

The symmetric algebra S(g) and the universal enveloping algebra U(g) are related 
through a map A, called symmetrization: 

A : S(fl)— W( fl ). (A4) 

Given that X and Y are elements of the Lie algebra q, XY will denote their product in 
S(g) (XY = YX), and X ■ Y will denote their product in U(g) (X ■ Y = Y ■ X + [X, Y}). If 
{Zi, Z 2 , . . . , Z n } is a basis for g, the symmetrization mapping is, by definition, 

A(Z?Z? . . . Z?) = ± Yl ■ Z ^ ■ ■ ■■ Z ^ » ( A5 ) 

where m = J2i e ^ *2, • • • , «m) is a set of m integers such that exactly €j of them are equal 
to j (1 < j < n), and the sum is over all permutations p of the m integers (1,2, . . . , m). 
For example, A(Z 1 Z 2 ) = \ (Z 1 ■ Z 2 + Z 2 ■ Zi). The symmetrization A is a vector space 
isomorphism (a canonical linear bijection) [?]. 

Since symmetrization is an isomorphism, its inverse is defined. Hence the extension of 
the domain of tp from the Lie algebra to the enveloping algebra is given by 

V^oA" 1 :U{g) ^C°°( *,M). (A6) 
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FIGURES 



FIG. 1. The deformation f3 is a unique function of the total angular momentum /. At the 
maximum angular momentum I = \ + fi, the deformation is a minimum [j m \ n = |A — /u|/3. 

FIG. 2. The triaxiality parameter 7 as a function of the total angular momentum / for rotation 
around the long and short principal axes when (A, fi) = (8,4). 

FIG. 3. For (A, fi) = (8, 4) the projection K onto the symmetry axis of the total angular 
momentum / of a rotating prolate spheroid is plotted versus the angular momentum. 

FIG. 4. For (A, fi) = (8, 4) the projection of the angular momentum on the short axis is plotted 
versus the total angular momentum for a sequence of triaxial tilted rotors which are critical points 
of energy functional I. 

FIG. 5. For = (8,4), sequences of triaxial principal plane rotors are shown as plots of 7 

and (Ii/I) 2 versus the angular momentum / when v = 20/3 for energy functional II. 

FIG. 6. For (A,/i) = (8,4), sequences of triaxial principal plane rotors are shown as plots of 7 
and (Ii/I) 2 versus the angular momentum / when v = for energy functional II. 

FIG. 7. For (A, fi) = (8,4), sequences of triaxial principal plane rotors are shown as plots of 7 
and (1 1 /I) 2 versus the angular momentum / when v = —4/3 for energy functional II. 
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TABLES 





TABLE I. 


Deformations for rotation around the 1-axis, fi = 




//A 


9i/A 


<? 2 /A 


93/A 


7° 



1 

4 
1 
2 
3 
4 
1 


l 

3 
1 
3 

_ 1 
3 
1 

3 
1 
3 


l 

3 

,1 VTE 
^6 8 

,i 

"•"6 4 
, 1 V7 

~6 8 

+ 1 


^3 

^6^8 
+ 1 + ^3 


0.00 
81 
3.69 
11.12 
60.00 




TABLE II. 


Deformations for rotation around the 1-axis, A = 




i/n 


?i/A» 




<?3/V 


7° 




1 

4 
1 
2 
3 
4 

1 


+§ 


+ \ 
1 , VTE 

6 "•" 8 
1 , V3 
6 4 
1 , %/7 
6 8 
1 
6 


2 
3 

1 %/i5 
6 8 

1 v 7 ! 

6 4 
1 77 
6 8 
1 
6 


60.00 
59.19 
56.31 
48.88 
0.00 



TABLE III. Tilted rotation of spheroids for the (8,4) su(3) representation. The projection of 
the angular momentum on the symmetry axis is K. 



Prolate Oblate 

1_ q K q K_ 

4 -3.33 4.00 

5 -3.22 4.04 

6 -3.07 4.07 

7 -2.89 4.10 

8 -2.67 4.07 2.67 8.0 

9 -2.39 3.96 2.39 8.3 

10 -2.03 3.62 2.03 8.7 

11 -1.54 2.34 1.54 9.4 

12 -.- -.- 0.67 12.0 
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m A FIT 1 — ' tt j 

TABLE IV. 


Principal plane rotations 


h = 


0, for the ( 


3,4) su(3) representation 




I 


?l 




% 


h 


9l 




% 


h 


4.00 


-3.33 


-3.33 


6.67 


4.00 


-.- 






-.- 


5.00 


-3.33 


-3.10 


6.44 


4.05 


-.- 


__ 


__ 


-.- 


6.00 


-3.33 


-2.80 


6.14 


4.11 


— 






— .— 


7.00 


-3.34 


-2.42 


5.76 


4.19 










8.00 


-3.35 


-1.92 


5.27 


4.31 


-5.33 


2.67 


2.67 


0.00 


9.00 


-3.37 


-1.24 


4.61 


4.48 


-4.76 


2.07 


2.68 


3.47 


10.00 


-3.54 


0.06 


3.48 


4.80 


-3.88 


0.93 


2.95 


4.78 


10.08 


-3.68 


0.50 


3.18 


4.84 


-3.68 


0.50 


3.18 


4.84 



TABLE V. X 3 


values for the ll K 


= 0" band of ( 


8,4). The values for q\,q2 


(?3 are calculated 


from Eq. (|3.18)) 












/ qi 


qi 


qz 


7° 




^3 










density 


eigenvalue 


-5.33 


6.67 


-1.33 


19.11 


0.00 


0.00 


2 -5.33 


6.54 


-1.21 


20.02 


-21.33 


-24.52 


4 -5.33 


6.13 


-0.80 


23.13 


-85.33 


-93.90 


6 -5.33 


5.31 


0.02 


30.20 


-191.98 


-214.23 


8 -5.33 


2.67 


2.67 


60.00 


-341.31 


-383.81 





TABLE VI. X 3 values for the "K 


= 4" band of prolate spheroids for 1 


8,4) 


I 


q 


density 


i"C3) 

eigenvalue 


4 


-3.33 


106.56 


108.18 


5 


-3.22 


77.17 


72.06 


6 


-3.07 


42.04 


31.94 


7 


-2.89 


4.13 


-16.23 


8 


-2.67 


-38.20 


-54.89 


9 


-2.39 


-81.20 


-123.64 


10 


-2.03 


-123.20 


-141.88 


11 


-1.54 


-161.04 


-250.00 
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TABLE VII. X 3 values for the U K = 2" band of (8,4). The fundamental system is solved for 
rotations in the 1-3 plane 



I 


Ql 


( l2 


Q3 


7° 


v(x 3 ) 

density 


eigenvalue 


2 


-5.06 


-1.60 


6.67 


16.70 


26.67 


24.52 


3 


-5.07 


-1.44 


6.51 


17.83 


0.69 


0.00 


4 


-5.08 


-1.21 


6.28 


19.57 


-35.78 


-14.28 


5 


-5.08 


-0.88 


5.96 


22.14 


-82.90 


-72.00 


6 


-5.09 


-0.43 


5.53 


25.96 


-140.86 


-70.71 


7 


-5.11 


0.21 


4.90 


32.04 


-210.35 


-183.10 


8 


-5.13 


1.32 


3.82 


44.27 


-292.52 


-159.29 


9 


g C 


G C 


G C 






-336.36 


10 


G C 


G C 


G C 






-291.45 



TABLE VIII. Deformations for Ne calculated with the cranked anisotropic oscillator and 
the density matrix method 



I Qi Q2 <?3 





cranking 


density 


cranking 


density 


cranking 


density 





-2.67 


-2.67 


-2.67 


-2.67 


5.33 


5.33 


2 


-2.67 


-2.67 


-2.58 


-2.54 


5.25 


5.21 


4 


-2.69 


-2.67 


-2.31 


-2.13 


5.00 


4.80 


6 


-2.72 


-2.67 


-1.78 


-1.31 


4.50 


3.98 


8 


-2.67 


-2.67 


1.33 


1.33 


1.33 


1.33 





TABLE IX. 


Deformations for 


the U K = 0" 


band of 24 Mg calculated with the 


cranked 


an: 


isotropic oscillator and Eq. (|3.18|) from the density matrix method 






I 




qi 




Q2 








cranking 


; density 


cranking 


density 


cranking 


density 





-5.33 


-5.33 


-1.33 


-1.33 


6.67 


6.67 


2 


-5.34 


-5.33 


-1.23 


-1.21 


6.56 


6.54 


4 


-5.34 


-5.33 


-0.90 


-0.80 


6.25 


6.13 


6 


-5.36 


-5.33 


-0.26 


0.02 


5.62 


5.31 


8 


-5.33 


-5.33 


2.67 


2.67 


2.67 


2.67 
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TABLE X. Deformations for the band heads of 24 Mg calculated with the cranked anisotropic 
oscillator and Eq. ( |3.19| ) from the density matrix method 



I 


qi 


<?2 




<73 






cranking density 


cranking 


density 


cranking 


density 





6.67 6.67 


-5.33 


-5.33 


-1.33 


-1.33 


2 


6.66 6.67 


-5.09 


-5.07 


-1.58 


-1.60 


4 


6.67 6.67 


-3.33 


-3.33 


-3.33 


-3.33 
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